Finite Temperature Magnetic Properties of Rare Earth Iron Garnets 
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Magnetic garnets as insulators are very important materials which have been studied We see that, as the supercell is cooled from 800K to approaching OK, it first 25 S We begin by applying DFT using the Vienna ab initio simulation package (VASP) to get 


extensively since the mid-twentieth century with promising outlooks in applications reaches the Curie temperature around 600K. This is the point at which the the magnetic parameters of the RIGs. We first write the spin Hamiltonian between two 


over many different areas of research. More specifically, Rare Earth - Iron Garnets moments of the tetrahedral and octahedral iron ions begin to be ordered, providing = 
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(RIGs) of the formula R3Fe5012 (Fig. 1) are -e a finite value magnetization of the system. Next, the compensation temperature 
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especially promising because of their low X 3 is reached at about 280K, which is the next point at which the net moment of 


Gilbert damping, high Curie temperature, IDA Vg e the supercell reaches zero. At this point, the gadolinium ions begin to be ordered in 
AD where the first term describes the collinear interactions between the ions, the second term 
and tunable magnetic properties [1]. They A AEN Son moments quite rapidly, leading to overtaking and quickly dominant a Temperature (K) 
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have been studied experimentally and Fih \en' influence of the net magnetic moment of the supercell. By the time the 
conversion [5]. 
theoretically in the past, but in order to temperature is approaching 0 K, the net moments of the gadolinium, octahedral iron, and tetrahedral irons are 21, 10, and -15 Bohr Magnetization 
Figure 1: R3Fe;O;2 structure: iron in brown, We can apply DFT as defined above on four ordered spin states to determine the J ij 
fully understand the magnetic properties of aa A per formula unit (upg/f.u.) respectively, leading to a net moment of 16 up/f.u. The total size of the supercell used was 4 x 4 x 4 unit cells of 160 
parameters [5]. Those states are as follows: 
these garnets, a computational atomistic approach should be applied. This will give, e.g., atoms each, totaling 10,240 atoms (6144 oxygen, 1024 octahedral iron, 1536 tetrahedral iron, and 1536 gadolinium.) 
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SO D The application of DFT to each of these states give us an energy value which we can 


It then becomes apparent that the application of the 


finally apply the following equation: 
RIGs have been found to have two different types 
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Which gives us the spin lattice coupling between any two ions in the supercell. It is an 


Figure 5: Magnetic moment of various elements in the supercell versus temperature for multiple values of applied magnetic field 


geometry (Fig. 2.) These different ions have been little effect on either of the irons. Both the 


important note to make that we use collinear magnetism to extract the magnetic 


found to have opposite magnetic moments of compensation temperature and Curie temperature were also found to have no significant change in response to a change in the magnitude of the l l l l l 
an parameters, implying that, under no applied field, there is no energy preference to any 
differing magnitude, meaning the garnets are applied field. It must be noted that the direction of the applied magnetic field is of no importance. No matter the direction of the applied field, the ae l 
Figure 2: Diagram of the bond angles of the direction of the total magnetic moment of the system. 
octahedral and tetrahedral irons [1] 
Jerrimagnetic. In addition to this, the rare-earth same results are found because of our adopted collinear magnetism. spene Heat apaciy Versus Temperatine l l 
13 NEPE From these electronic configurations, we re-scale the couplings and run various Monte- 
ions have been found to have a drastic increase in their magnetic moment below a We then look at the specific heat capacity versus temperature for 16 0395 we a l 
0 Oe io Carlo heat bath simulations on the data. Using the Heat Bath Monte-Carlo code [6], We 


specific temperature, called the compensation temperature, Ty. This moment is in the the system under different applied fields (Fig. 6). Note that the large l B 
can look at many different properties of the RIGs under many different conditions. 
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same direction as that of the octahedral iron, and, at low enough temperatures, stronger at 583K is associated with the phase transition at the Curie au o 
i Cs The instance of the heat bath code that we run allows us to look at the magnetic properties 
in magnitude than both iron ions [1]. Using the equation: temperature. As the applied magnetic field at the compensation oe ce ee . l Saa 
© 207 Applied Fed of the garnets under a variety of scenarios. We begin with the Monte-Carlo code by first 
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ee oC o4 running the same zero field cooling (ZFC) computations that have been done both 
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l n l l l l scaling factors until the calculations match the experimentally found results. From there, 
Figure 6: Specific Heat of the supercell versus temperature under different magnitudes of applied magnetic 
decreases to 0 K, one can find a compensation temperature at which the RIG has zero In order to ascertain that there is a phase transition, we next field , , , i ' . 
P P we look at the RIGs under field cooling. Field cooling (FC) is the experimental procedure 
net magnetic moment [2]. This temperature has been found to coincide with the analyze the magnetization vectors of the supercell across all temperatures and all applied magnetic fields. l an 
of applying a magnetic field to the compound while at the same time cooling the material 
relatively low temperature at which the rare-earth elements begin to have drastic Cosine of Angle from Z — Axis of Octahedral Iron Magnetization Vector Figure 7 plots the angle from the Z - axis of the octahedral iron 


down, as opposed to the ZFC mentioned before. We use this procedure to see the effect of 


changes in magnetic moment [1]. It must be also noted that the oxygen atoms in the magnetization vector as a function of temperature for different levels of 


the application of a magnetic field on the magnetic moments of each individual atom, and 


supercell have no contribution to the magnetic properties as far as the extent of this applied field. At OT applied field, the direction of the magnetization 
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paper is concerned. vector with respect to the Z - axis is seemingly random, as consistent 


The calculations were carried out on the University of Arkansas high performance 


with collinear magnetization. This is because of the lack of Zeeman : , , 
computing cluster, Razor. All the data was then processed using Microsoft Excel, and all 


effect - which states that the total magnetic moment should align with 


figures were made in either the modeling software Vesta [7] or in Microsoft Excel. 
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Because of the relatively large size of these RIGs (160 atoms per unit cell), few first- 
Figure 7: Cosine of the angle from the Z-axis of the unit vector in the direction of the net magnetization for 
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principle calculations were carried out on them to determine properties until reference 


direction. With an applied field, however, the Zeeman effect is demonstrated, and the total net magnetic moment is in the direction of the applied 


magnetic field ([001] direction). As can be seen in the graph, once the temperature hits the compensation temperature, the magnetization of the 
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Figure 3: [a] Experimentally found values for the net magnetic moment of various RIGs as a function of temperature ‘ ; | 
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In this study, we further investigate the magnetic properties of a rare-earth garnet to a decreased past the compensation temperature. The red arrows indicate the relative 
magnitude and direction of the moments of the ions. 
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the system, as schematized in Figure 8. 
better understanding of the magnetic properties of RIGs, and we possibly look to 


uncover any new phenomena that may be possible with them. 
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